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SATURATION  OF  THE  LOWER-HYBRID-DRIFT  INSTABILITY 
BY  MODE  COUPLING 


A  first  principles  theory  of  the  nonlinear  saturation  of  the  lower- 
hybrid-drift  instability  and  associated  particle  transport  is  presented.  A 
nonlinear  equation  describing  2-D  lower-hybrid-drift  turbulence  is  derived 
which  includes  the  nonlinear  E  x  B  and  polarization  drifts  of  the 
electrons,  adiabatic  ions,  and  realistic  sources  and  sinks  of  energy.  The 
equation  is  solved  numerically  to  obtain  the  evolution  and  saturation  of 
the  wave  spectrum,  as  well  as  the  self-consistent  particle  flux  which 
yields  the  diffusion  coefficient. 

The  lower-hybrid-drift  instability  is  a  flute  mode  (k  •  B  ■  0)  which 
is  driven  unstable  in  a  plasma  with  strong  inhomogeneities,  p./L  > 

«  /A  “  ® 

(m  /m. )  ,  where  p  and  L  are  the  ion  Larmor  radius  and  density  scale 

e  i  in 

length,  respectively.  The  linear  behavior  of  the  Instability  is  well 

understood. In  the  weak  drift  regime  (V^  <  v^,  where  *  vi^niLn  *8 

the  ion  diamagnetic  drift  velocity  and  v^  is  the  ion  thermal  velocity)  the 

mode  is  driven  unstable  by  the  resonant  interaction  of  the  wave  with  the 

drifting  ions  when  u>/k  <  V  .  The  growth  rate  maximizes  at  k  -  p  1  = 

-i  1/2  y  <*1  7  es 

p .  (m,  /m  )  «  The  mode  is  stable  for  k  p  »  1,  and  approaches  marginal 

i  i6  ~  y 

stability  as  k  p  ♦  0.  The  wave  frequency  is  such  that  0.  «  m  «  0  so 
y  68  X  € 

that  the  electrons  are  strongly  magnetized  while  the  ions  can  be  treated  as 
unmagnetized.  In  a  finite  B  plasma  the  resonant  interaction  of  VB  drifting 
electrons  with  the  wave  is  stabilizing. 

The  nonlinear  behavior  of  the  instability  is  much  less  well 
understood.  Ion  trapping  has  been  observed  to  quench  the  growth  of  the 

4 

instability  in  particle  simulations.  The  onset  of  stochastic  electron 
heating  has  also  been  proposed  to  explain  saturation  in  these 
simulations.^  On  the  other  hand,  ion  trapping  is  not  a  viable  saturation 

mechanism  when  a  broad  2-D  spectrum  of  waves  is  excited,  and  it  has  not 

* 

been  observed  in  recent  simulations  using  realistic  mass  ratios  where  such 
spectra  develop.^  In  this  letter  we  focus  on  the  nonlinear  coupling  of 
stable  and  unstable  waves  to  saturate  the  Instability. 


Ue  consider  a  slab  equilibrium  of  cold  electrons  and  warm  ions  with  a 
density  profile  Oq(x)  supported  by  a  magnetic  field  Bzq(x).  The 
equilibrium  ion  velocity  distribution  is  taken  to  be  a  Maxwellian  with  an 

A 

average  drift  velocity  e^.  Because  of  the  flute  nature  of  the 

instability  (k  •  B  ■  0),  we  self-consistently  limit  the  spatial  variations 
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to  the  xy  plane  (as  contrasted  with  universal  drift-wave  turbulence  which 
Is  Inherently  3-D).  In  the  weak  drift  regime  the  ions  respond  to  the 
perturbed  potential  <Kx,y,t)  adlabatically  to  lowest  order  since  to  «  kv^. 
Thus, 

ni  “  no^ 1  +  /7  Vl1,v”1|^lr  +  vdi  l7)]exP<-e^Ti>  <l> 

where  the  term  proportional  to  /tT  in  (1)  is  a  small  correction  describing 
the  resonant  ion  interaction.  The  perturbed  electron  motion  is  simply 
given  by  the  E  x  B  and  polarization  drifts, 


V 

~e 


x  e 

z 


(2) 


where  d/dt  »  3/3t  +  V  »7. 

~e 

be  written  as 


Using  (2),  the  electron  continuity  equation  can 


(3) 


Finally,  invoking  charge  neutrality  (ne  -  nt)  and  combining  (1)  and  (3),  we 
obtain  the  nonlinear  equation 

a  a  a  — «  a  a  a 

(i  -  7  )$  +  <p  -  yJv  l(d  +  *  )  +  Y  4> 

T  y  0  ytt  rty  eY 


+  (7$  x  e  »7)  7%  +  y  v}  x  e  •  7|7'1|(J  +  l)  -  0 

z  o  z  x  y 


(4) 


where  7  -  p  7,  t  -  (p  /L  )Q  t,  $  -  (e$/T  )(L  /p  )  and  the 
.  es  es  n  e  1  n  es 


subscripts 

on  $  denote  a  derivative  with  respect  to  that  variable.  The  quantity 

represents  the  wave  damping  due  to  7B  resonant  electrons  in  a 

1/2 

finite  8  plasma.  Equation  (4)  is  only  valid  for  ya  ■  »  7.,/v.  <  1  since 
the  adiabatic  ion  response  can  only  be  justified  in  this  limit. 
Linearizing  this  equation,  we  obtain  the  complex  eigenvalue  (in  our 
normalized  units), 

A  AAA 

-  -  K  k‘|k|  . 

-i  -  iY  — rr  . 

2x3  e  ^ 2 


-  0), 


*  j  +  iY  n 

1+k*  ( l+k“) 


u 


(5) 


In  Che  Holt  Y  ,  yn  *  0,  (4)  reduces  to  the  Hasegawa-Mima  equation  in 
which  the  nonlinearity  arises  from  the  nonlinear  polarization  drift.  This 
equation  has  two  invariants,  energy  and  enstrophy,  neither  of  which  is 
preserved  in  the  more  general  (4).  When  Yq  is  finite,  the  E  x  B 

nonlinearity  also  appears  in  (4).  This  nonlinearity  has  been  considered  by 

a  q 

Horton  and,  more  recently  by  Waltz,  in  studying  universal  mode 
turbulence.  However,  our  equation  differs  from  their  work  in  that  we  do 

A  ^  A  A  A 

not  assume  that  $  +  <p  is  replaced  by  it  linear  counterpart  -i(o>  -  k  )$. 

t  y  k  y 

Thus,  our  equation  is  second  order  in  time  rather  than  first  order. 

The  direction  of  energy  flow  described  by  (4)  can  be  understood  by 

A 

calculating  the  stability  of  a  single  large  amplitude  wave  with  k^ 

and  u>q  satisfying  (5).  A  perturbation  of  wavevector  k  is  coupled  through 

the  pump  to  modes  with  k  +  pk.  (p  ■  1,2,...).  For  simplicity,  we  consider 

«  **  —  * 

only  the  coupling  of  (u>,k)  with  its  nearest  neighbors  (m  ,k  )  where 

A  A  A 

o»+  -  u>  +  o)q  and  k.^  ■  ]c  +  kg,  !•«*,  P"l.  The  dispersion  relation  for 
this  (k.Jt^)  coupled  system  is 

M(k  ,L)M(k,L) 

+  H0\Z[ - v  — -  +  <fe0  ♦  -k^0)]  -  0  (6) 

where  e(k,u>)  -  <o(l  +  k2(l  -  id  )J  -  k  +  iY  ,  6.  -  yQ(kv-u)/k3  and 

y  a  “  k 


Y0(ky““)/k  ^ 


M*ki,k2^  "  -1  X  V®z  Ikl(1  “  i6kl)  "  k2jl  "  \6k2)h  WheD  Y0  and  Ye  are 

neglected  in  (6)  and  we  take  the  limit  u  »  w  ,  the  decay  modes  are  purely 

a  A  A  A 

growing  with  a  growth  rate  which  peaks  around  k  •  kn  *  0  with  |k|  <  Ik  | . 

A#  -WQ  A#  A*Q 

A  necessary  requirement  for  instability  is  that  one  of  the  decay  waves 

A  A 

(k,k|.)  has  a  longer  wavelength  than  the  pump.  For  this  situation,  in  which 
(4)  reduces  to  the  Hasegawa-Mima  equation,  the  wave  energy  inevitably 
cascades  to  longer  and  longer  wavelength  so  that  no  stationary  wave 
spectrum  can  result. When  Yn  is  included,  this  conclusion  no  longer 

»  A  A  A  A 

remains  valid.  Taking  the  limit  Ik^l  «  |k|  and  again  assuming  u  »  o>^, 
the  dispersion  relation  is  given  by 

A  A  A  A  A  A  A  A  A  A 

*l*0l  k„  .  .  V<2*  >, 

ur  ■  *2  L  ^  *  3  "2  '* 

•  i  +  k  Ikr(i-Ht^) 

Equation  (7)  yields  a  dissipative  instability  which  is  driven  unstable  by 
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the  E  x  B  nonlinearity.  This  instability  produces  a  flow  of  energy  from 
long  to  short  wavelengths  and,  as  will  be  demonstrated,  enables  the  wave 
spectrum  to  reach  a  steady  state. 

We  solve  (4)  numerically  by  decomposing  it  into  two  coupled,  first- 
order  (in  time)  differential  equations  which  are  advanced  in  time.  These 
equations  are  solved  by  a  spectral  method  developed  by  Fyfe  et  al.,^  based 
on  the  work  of  Orzag.^2  The  potential  $  is  Fourier-decomposed,  i.e., 

*  A  A 

♦  ~  exp(ik  •  3)  where  k  -  n/X  and  g  «  n  e  +  n  e  with  n*  and  ny 

x  x  y  y  / 

integers,  and  X  -  5  (which  fixes  the  |n|  for  which  |k|  -  1).  The  numerical 
results  presented  in  this  letter  are  nominally  computed  on  a  32  x  32  mesh. 

The  electron  damping  y^  controls  the  region  of  instability  in  k  space 
by  stabilizing  short  wavelength  modes.  We  have  chosen 

Ye  •  Y0!k|7(1  +  k2)"V4  (8) 

with  •  1.6,  which  yields  an  unstable  wave  spectrum  which  is  realistic 
for  the  lower-hybrid-drift  Instability  in  finite  8  plasmas.  The  growth 

AM  A  A 

rate  peaks  at  (k  ,k  )  *  (0,+l)  and  is  stable  for  k  >  1.6  and  k  >  1. 

x  y  —  y  x 

A  A  ^2  — 2 

The  spectrum  is  initialized  with  random  noise  with  $  -  10  -  10  , 

and  the  system  is  allowed  to  evolve  until  the  wave  energy,  given  by 

W/nT  -  (me/xm1)  yJ  E  (1  +  k2)  lij2  (9) 

k 

reaches  a  steady  state  value.  The  mean  square  of  the  potential 
fluctuations, 

P  -  <(e*/T)2>1/2  -  [(2m  /wm)  yl  Z  ij  \2]lJ2  (10) 

e  k 

as  well  as  the  nonlinear  and  quasilinear  diffusion  coefficients 

A 

*  1/22*'*  **  *  * 

D  -  D  ./v  p  -  (2m  /wm.)  '  yf  Z  (k  /|k|)  <f>.  (k  <p,  -  i~)  (11) 

nl  nl  es  as  e  i  0  ^  y  *<  7  Tk  '  y  k  3t  7  v  7 

*  1/9  9  *  *9  *  9  *9—1 

°,1  ■  VV..PM  ■  ’'O  l  k,  ^k1  +  k )  <l2> 

respectively,  are  also  computed.  The  nonlinear  diffusion  coefficient  d 
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Is  derived  by  averaging  the  particle  flux  in  the  x  direction  over  y,  i.e., 

A 

«  <nV^>^.  The  quasilinear  expression  is  then  obtained  by  approximating 
3/3t  -  -im^  an<*  u8ln8  (5)  in  (11). 

In  Fig.  1  we  show  Instantaneous  2-D  wave  spectra  in  n  space  for 
Yq  ■  .5.  Figure  la  shows  the  spectrum  during  the  linear  phase  of  the 
instability,  and  Fig.  lb  shows  the  spectrum  after  saturation.  In  the 
linear  regime,  the  wave  spectrum  is  strongly  peaked  around  the  most 
unstable  modes,  (n^ny)  ■  (0,+5).  The  spectrum  begins  to  saturate  as  these 
modes  couple  through  a  n  ■  0  mode  (typically  with  n  ~  4)  to  damped  modes 

y  * 

with  nx  £  ny.  A  peak  in  the  wave  spectrum  appears  around  (4,0)  at  this 
time.  This  transient  phase  is  completed  as  the  total  wave  energy  saturates 
and  spreads  through  most  of  the  unstable,  weakly  damped  volume  of  n  space 
as  shown  in  Fig.  lb.  The  spectrum  typically  remains  peaked  around  nx  -  0 
with  ny  typically  somewhat  smaller  than  that  of  the  linearly  most  unstable 
mode  [see  Fig.  lb  for  which  the  dominant  modes  are  (n  ,n  )  -  (0,+3)].  The 
2-D  wave  spectrum  exhibits  substantial  variability  in  time,  even  after 
saturation,  as  the  unstable  and  stable  modes  continue  to  exchange  energy  in 
a  dynamic  manner. 

The  time  evolution  of  the  total  wave  energy  (W),  mean  square  of  the 

A 

potential  fluctuations  (F),  and  the  nonlinear  (D  ,)  and  quasilinear 
a  ni 

diffusion  (D^)  coefficients  are  shown  in  Fig.  2.  All  of  these  quantities 
exhibit  a  similar  temporal  behavior.  The  initial  decay  (t  <  4)  is 
associated  with  the  rapid  dissipation  of  energy  initialized  in  the  damped 
modes,  and  is  followed  by  a  linear  growth  phase  (4  <  t  <  20).  Subsequently, 
mode  coupling  occurs  which  leads  to  saturation  of  the  instability,  albeit 
with  some  initial  overshoot  (t  ~  24).  The  levels  of  the  total  wave  energy 
and  other  parameters  of  Fig.  2  are  quite  stationary  in  time  after 
saturation.  Also,  the  stationary  values  of  all  four  quantities  are 

A 

relatively  insensitive  to  the  initialization  of 

A 

We  note  that  in  Fig.  2,  the  quasilinear  diffusion  coefficient  (D  , ) 

A 

tracts  the  actual,  nonlinear  diffusion  coefficient  (D^)  quite  well  during 
the  entire  time  evolution  of  the  instability.  An  important  point  which 
must  be  emphasized  with  regard  to  is  that  both  species,  electrons  and 
ions,  continue  to  exchange  both  energy  and  momentum  even  after  a  steady 
state  is  reached;  the  electrons  through  the  VB  resonance  and  the  ions  by 


direct  resonant  Interaction.  If  the  Instability  had  saturated  by  Ion 
trapping  and  the  electrons  had  no  resonant  Interaction  with  the  wave,  there 
would  be  no  diffusion  In  the  steady  state  since  the  electrons  could  not 
exchange  momentum  with  the  Ions.  Both  species  oust  have  a  dissipative 
interaction  with  the  waves  to  have  diffusion. 

A  number  of  runs  have  been  made  for  different  values  of  ‘v  ,  the  drift 

.  0 

parameter.  The  saturation  level  of  $  is  found  to  be  nearly  Independent  of 

*  2  2  2 

Yq*  Thus,  from  Eqs.  (11)  and  (12),  It  is  found  that  «  Yq  *  V^/v^. 

* 

Next,  we  compare  our  stationary  values  of  P  and  D  ,  with  these 

v»  nl 

obtained  from  “wave  energy  bound"  calculations.  The  "energy  bound" 

results  from  equating  the  wave  energy  with  the  kinetic  energy  associated 

2 

with  the  relative  drift  of  the  electrons  and  ions,  men(V^£  -  Vj#)  /2.  This 

free  energy  bound  Is  actually  an  extreme  underestimate  of  the  free  energy 

in  a  finite  6  plasma,  since  the  instability  can  also  feed  off  the  magnetic 

free  energy,  which  typically  greatly  exceeds  the  drift  energy. ** 

Nevertheless,  this  bound  has  been  widely  Invoked  and  seems  to  yield 

reasonable  transport  rates  when  compared  with  experimental  observations. 

1/2 

The  "drift  energy  bound"  yields  a  potential  Pab  *  (m^/wm^  yq,  where  we 
have  assumed  Te  «  T^.  The  scaling  of  Peb  is  identical  with  that  given  In 
(10)  since  $  is  Independent  of  Yq.  For  the  case  Yq  *  .5,  P/Pej,  *  3  (from 
Fig.  2)  so  that  the  actual  fluctuation  level  obtained  from  our  code  is 
slightly  larger  than  calculated  from  the  "energy  bound".  Finally  we 

A  A 

compare  D  ,  with  that  obtained  from  the  simple  formula  D  -  (Y/k2)  /V  p  , 
r  nl  r  'm,  es  es’ 

where  Y  and  k  are  evaluated  where  Y  peaks.  We  find  D"  (m  /irm. )  y«. 

1  e  1  0* 

which  again  has  the  same  scaling  as  Dq^  but  is  a  factor  of  2.4  smaller  for 
Yq  *  0*5  (from  Fig.  2). 

Measurements  of  lower-hybrid-drift  turbulence  in  a  9-plnch  by  CO2 
laser  scattering  have  been  reported  recently. The  observed  wave  spectra 

A 

were  flute-like  (k  •  B  ■  0)  with  clear  peaks  around  k  *  0.7  (a  factor  of 
^  ^  y  .  . 

two  smaller  than  the  most  unstable  linear  mode),  and  k  <  k  .  Data  was 

x  y 

taken  for  three  different  filling  pressures,  corresponding  to  three  values 

of  V41/V1.  Electron-ion  collisions  were  apparently  significant 
1/2 

(v  ^  -  (fl  0^)  ]  for  the  two  highest  filling  pressures  so  we  can  only 

compare  our  results  with  the  data  from  the  lowest  filling  pressure,  which 
corresponds  to  V^j/v^  ■  0.52  (Yq  3  0.9).  The  measured  density  fluctuation, 


n/tig  -  0.014,  is  quite  close  (smaller  by  a  factor  of  2)  to  our  theoretical 
prediction  of  n/n^  3  0.025.  In  addition,  the  shift  to  long  wavelength  is 
consistent  with  our  calculated  spectra  (compare  Figs,  la  and  lb). 

Computer  simulations  of  the  lower-hybrid-drift  instability  with 
realistic  mass  ratios  have  also  been  performed  recently  using  a  fluid- 
kinetic  numerical  scheme.^  Unfortunately,  these  computations  have  been 
carried  out  for  V^/vj^  >  1,  which  is  outside  the  range  of  validity  of  our 
present  theory.  Nevertheless,  a  broad  spectrum  of  modes  is  observed  in  k 
space  in  the  simulations  which  is  consistent  with  our  calculations. 

In  conclusion,  the  lower-hybrid-drift  instability  can  saturate  via  a 

mode  coupling  process.  The  basic  saturation  mechanism  involves  a  transfer 

of  energy  from  the  growing,  long  wavelength  modes  to  the  damped,  short 

wavelength  modes.  The  nonlinear  E  x  B  drift  plays  a  crucial  role  in  this 

saturation  mechanism  by  preventing  the  polarization  drift  nonlinearity  from 

driving  the  wave  energy  towards  the  undamped,  long  wavelength  modes.  The 

saturation  energy  and  diffusion  coefficient  associated  with  the  wave 

2 

turbulence  scales  as  (V^/v^)  ■ 


2  2.6 
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